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Abstract. The evolution of central stars of planetary nebulae can proceed in several distinct 
ways, either leading to H-deficiency or to H-normal surface composition. Several new simulations 
of the evolution channels that lead to H-deficiency are now available, mainly the born-again sce- 
narios that are triggered by a He-shell flash during the hot pre-white dwarf evolution phase. 
A realistic AGB progenitor evolution is important for correct HRD tracks, that allow mass de- 
terminations. New hydrodynamic simulations of He-shell flash convection including cases with 
H-ingestion are now performed, and allow a determination of the convective extra-mixing effi- 
ciency. This has direct consequences for the intershell abundance distribution of AGB stars that 
can be observed in the H-deficient CSPN. 
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1. Introduction 

The evolution of the central stars of planetary nebulae (CSPN) is closely connected 
to the evolution of the progenitor phase, the Asymptotic Giant Branch (AGB) stars, 
and has observable effects for the CSPN progeny, the WD stars. The AGB evolution 
phase is the result of previous main-sequence and horizontal-branch evolution of low- 
and intermediate mass stars. An up dated c l assific ation of mass ranges involving AGB 
stars has been recently suggested bv iHerwiel l|2005[ Fig. 2). Accordingly, low-mass stars 
are those below an initial mass of « 1.8 that ignite a He-core flash at the tip of 
the RGB, while intermediate mass stars are those initially more massive than « 1.8 Mq 
but less massive than required to explode as a core-collapse supernova. The dividing 
mass at « 10 Mq separates those super- AGB st ars that form ONe- white dwarfs and 
those that ignite as core collapse supernova (e.g. ICarcia-Berro fc Ibenlll994l) . However, 
quantitatively the fraction of AGB stars that explode as core-collapse supernova is still 
very uncertain, largely due to uncertainties in super- AGB mass loss, and the physics of 
mixing and burning, e.g. during dredge-up and hot-bottom burning (see below). For AGB 
stars it is useful to use the limiting mass for hot-bottom burning (HBB) as a dividing line 
for low-mass and massive AGB stars. Hot-bottom burning t ransforms dredged up C into 
N an d prevents the formation of C-rich composition fe.g. IScalo et al.lll97!Tl: iLattanziol 
Il992l) . HBB has therefore an important influence on the evolution of AGB star, and 
the composition of the ejecta that later form the planetary nebula. AGB stars that are 
massive enough to ignite carbon burning are the super- AGB stars. 
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Figure 1. Comparison of three different evolution channels in the g-T^ff-plane for a 0.604 M© 
post-AGB star that had a main-sequence mass of 2 M© and evolved through about a dozen 
thermal pulses with dredge-up. This figure originally appeared in the Publications of the Astro- 
nomical Society of the Pacific (Werne r fc Herwie 20 0fi'l. Copyright 2005, Astronomical Society 
of the Pacific; reproduced with permission of the Editors. 



Conceptually, the evolution of AGB stars is well understood, at least at solar and mod- 
est metal deficiency, like in the Magellanic Clouds. The review bv llben fc Renzinil l)l983f) 
describes the basic properties of AGB stars: thermal pulses, the necessity of the third 
dredge-up, the ef fect of mass loss and the poten tial f or impor t ant nu cleosynthesis. More 
recent reviews bv lLattanzio fc Boothrovdl l(l997^ and lBlockeil l|l99?t l include accounts of 
hot-bottom burning, and an improved quantitative description of AGB evolution. D edi- 
cated review s were written on mass loss (jWillson ,2000,). the s process (Busso et al. 19991: 
lMevejll994) a nd the cool post-A GB stars that allow important observations concerning 
the s process l|van Wincke]ll20f)l . 

Unfortunately, major unreso lved issues rem ain, that severely limit the predictive qual- 
ity of current AGB modeling ^|Herwigl|200l^^ . The issues concern the mass-loss history, 
the nature and importance of hydrodynamic mixing, including the physics of mixing 
at the bottom of the envelope convection zone that seems to be critical for the third 
dredge-up, and for formation of the assumed ^^C-pocket. Another problem is the rather 
poorly known effect of rotation and magnetic fields on the evolution, nucleosynthesis 
and formation of asymmetries of bi-polar structure in proto-planetary nebulae. All these 
problems are amplified for stars at very low metallicity, partly because less observational 
checks are available. 

The evolution of CSPN (described in Sect. 13) depends sensitively on the progenitor 
AGB evolution. In particular, the observed abundances of the H-deficient CSPN are 
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linked to mixing processes in the He-shell flash of the AGB phase, and some new results 
in this area are presented in Sect.|3| 



2. Evolution of CSPN 

We have now modern stellar evolution m odels of sev eral distinct channels of post- AGB 
evolution, and through the CSPN phase ijWerner fc Herwia 12006) . These channels are 
characterized by the thermal pulse cycle phase at which they depart from the AGB. The 
most common is the undisturbed evolution off the AGB, through the CSPN phase and 
onto the cooling track of white dwarfs. These CSPN should show a surface abundance 
pattern that reflect s the former AGB envelope including enrichment of that envelope 
through dredge-up I Napiwotzkilll999|) . 



T he possibili t y of a born-again evolution was recognized early, bv ISchonberned lll983l) 
andllben et alJ lll983l). Several ni ore recent calculations are now available ijHerwig et al.l 
ll999HLawlor & M acDonald 2003), inclu ding a new grid of H -deficient CSPN tracks based 
on the born-again scenario by Mille r Bertolami et al.l l)2006^) . Two flavours of the born- 
again evolution are now distinguished: in the very late thermal pulse (VLTP) the re- 
maining small amount of H in the surface layers is convectively mixed into the He- shell 
flash convection, wh ile in the late thermal pulse (LTP) this does not happen ((B locked 
l200lUHerwi^l2001bl) . It should be noted that the physics of convectiv e-reactiy e mixin g 
which characterizes the VLTP is not well understood, as emphasized bv lHerwidl|2001ah . 
Currently, efforts are underway at Los Alamos and elsewhere to study the hydrodynamic 
aspects of this problem. The simulations discussed in the following section are a first step 
into that direction. 

According to their different nucleosynthetic and mixing evolution the VLTP and LTP 
cases will produce CSPN with observationally distinguishable properties. For example, 
while Sakurai's object and V605 Aq are now considered the result of a VLTP evolution, 
PC Sge is the result of a LTP evolution, and t he arguments in the liter ature leading to 
this conclusion have been reviewed recently bv lWerner fc Herwid l|2006() . 

In addition to the two born-again channels there may be a possibility for a H-deficient, 
but not H-free CSPN evolution without the born-again sequence of events. The AGB 
final thermal pulse (AFTP) scenario requires some mechanism to shed the remaining 
small envelope mass exactly after the final AGB thermal pulse with efficient dredge-up. 
In view of the ph ysics uncertaintie s , in par ticul ar of envelope convectio n, this scenario is 
not inconceivable. IPe Marco et al.l 1 2008hl) and lPe Marco et alJ 1 2008al) have tentatively 



investigated the spiraling-in and swallowing of a massive planet or very low mass stel- 
lar companion as an enhanced mass-loss process for th e AFTP. The AFTP sc enario is 
currently the most likely scenario for the hybrid CSPN l)Napiwotzki et al.lll99l[) . 

Fig-fflshows the comparison of different evolution channels. These calculations give an 
estimate of the difference of the tracks due to different surface chemistry and evolution 
history. In particular the highest temperature is the same for all channels. 

More important for the location of the track in the HRD, and therefore for mass 
determinations using stellar-evolution tracks, is the progenitor evolution. Fig.|21 shows 
CSPN from various authors. The more recent 0.604 M© track is hotter than older tracks 
of the same mass. This difference is largely due to differences in the assumptions a bout 
AGB mass loss and the treatment of third dredge-up. As already shown by|B16cker (19951 ) 
CSPN with the same core mass can have significantly different tracks in the HRD. A core 
descending from an initially more massive progenitor simulated with a higher mass-loss 
rate will follow a hotter track than a core of the same mass with a lower-mass progenitor 
with lower mass loss. Similarly, efficient third dredge-up leads to a hotter CSPN track 
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Figure 2. H-nomial evolutionary traclts (labels: mass in M©) are from lSchon berne? ("1983?) and 
(Blocker 1995) (dashed lines) and Wood & Faulkner (1986). The 0.604 Mq track Herwie (20Q1); 
[Werner & Herwie (2006) (dot-dashed line) is the final CSPN track following a VLTP evolution 
and therefore has H-deficient composition. Shown as well are hot hydrogen-deficient post-AGB 
stars in gn^T^-plane. Wolf-Rayet central stars of early and late type ([WCE], [WCL], from 



iHamannl Jl997l^). PG1159 stars as well as two [WC]-PG1159 transition objects (Abell 30 and 



Werner fc Herwi jl2006l for ref.) Tee for the [WC] stars is related to the stellar radius 



78) (see 

at rRoss=20. This figure originally appeared in the Publications of the Astronomical Society 
of the Pacific (Werner & Herwie 2006). Copyright 2005, Astronomical Society of the Pacific; 
reproduced with permission of the Editors. 



compared with a core descending from a lower-mass progenitor with less third dredge-up. 
Binarity may influence the relation between core mass and HRD CSPN track additionally, 
through common envelope evolution, tidal synchronization, enhanced mass loss etc. 



3. Hydrodynamic simulations of He-shell flash convection 

We used the RAGE code to perform hydrodynamics simulations of a convection zone 
under conditions closely resembling those encountered in a typical He-flash during the hot 
pre-white dwarf evolution phase. The initial stratification is in hydrostatic equilibrium 
and consists of an entropy plateau, where the convection zone is expected to develop, 
with two adjacent stable regions. The convection is driven by a heat source in a thin 
layer at the bottom of the plateau to simulate the effect of He burning. An number of 2D 
and 3D runs with various resolutions and heating rates was produced. See lHerwig et all 
for more details. 

The first panel in Fig.|31shows the spatial entropy inhomogeneities at a snapshot during 
the onset of convection. Hot plumes rise from the bottom of the entropy plateau. They 
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Figure 3. Entropy inhomogeneities for three snapshots of a 2D simulation with a resolution 
of 2400x800 grid points. The horizontal average of the entropy has been subtracted to render 
small fluctuations visible. A brighter color indicates material with an entropy excess (i.e. lower 
density than the surroundings). The overturning flow inside the convection zone clearly differs 
from the oscillatory motions due to gravity waves inside the stable layers. 

merge while they grow in ampUtude and vertical extent and start to excite gravity waves 
in the upper stable layers well before actually reaching the upper boundary. Later, there 
are typically four large instationary convective cells wit h a number of smaller eddies that 
continue to drive gravity waves in the adjacent layers ijPres 
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Figure 4. Logarithm of the diffusion coefficient versus height for models with different resolu- 
tions (indicated in the legend) . The two vertical lines mark the boundaries of the entropy plateau 
at the start of the simulations. Later, they are slightly smeared out. The results from the runs 
with the two lowest resolutions differ significantly from the runs with next higher resolutions. 
Nevertheless, the agreement between the two highest-resolution runs is rather good. 

The transition between convectively unstable and stable layers appears to be very 
sharp - without any plumes penetrating significantly into the the stable region. This is 
very different from the situation of thin stellar surface convection zones (' Frevtag et al.l 
[l996 ). However, /c-cj-diagrams for the vertical velocity at various height levels clearly 
show the mode signature of the gravity waves well inside the convection zone an d also 
evidence of the convective flow reaching into the stable layers l|Herwig et al.ll2006|) . 

Of particular interest are the mixing properties across the convection boundaries. The 
mixing is efficient inside the convection zone. Sometimes blobs of relative hot material 
are sheared of the boundary to the top stable zone (see e.g. the tiny bright spots near 
the top left of the convection zone in the middle and bottom panels in Fig. and are 
transported rapidly into deeper layers once they are entrained in a downflow region. 
Future simulations that include various species and nuclear reactions will show what 
happens if H-rich material is brought in this way into layers where it can burn. 

The velocities in the stable layers are not much smaller than in the convection zone. 
However, wave motions are almost reversible and do not contribute much to the mixing. 
To quantify the mixing the trajectories of tracer particles (150 particles per horizontal 
layer in 68 layers) have been integrated based on snapshots with a high sampling rate of 
0.5 s. The spread of the relative vertical positions - or alternatively the relative entropy 
- of the particles in each ensemble with time gives a measure of the diffusion coefficient 
with height, plotted in Fig. 0] for models with different numerical resolutions. The large 
diffusion coefficient within the convection zone drops rapidly near the boundaries. This 
decay reaches slightly into the stable regions and can be described by one (at the top) 



The determination of the very small mixing coefficients far away from the convection 
zone is not very reliable. This small amount of mixing in the stable layer is due to internal 
gravity waves, and may be far reaching. 
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4. Conclusions 

Among the CSPN channels, the evolution leading to H-deficient surface composition is 
of great importance for constraining nucleosynthesis and mixing in the AGB progenitors. 
A clear correlation exists, e.g. between the mixing efficiency at the bottom of the He-shell 
flash convection and the O abundance in the intershell that is subsequently observed 
in H-deficient CSPN. AGB calculations with parameterized overshoot calculations by 
iHerwid l|2000l) showed that more convective extra-mixing leads to higher O abundance, 
and that an exponential decay parameter for extra-mixing (see Jlerwig (2000) for the 
definition of) / w 0.01 would be consistent with CSPN observed O abundances in the 
range 8 — 20% by mass. Our preliminary analysis of hydrodynamic simulations of the He- 
shell flash convection zone quantitatively agrees with this level of convective extra-mixing. 
We regard this agreement of hydrodynamic simulation and astrophysical observation as 
a successful validation of the computational simulations, and feel encouraged to continue 
this work towards more complex simulations of convective-reactive mixing. 
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